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Abstract 

The Minimal Prescription procedure is applied to tame the Landau pole singularities of 
resummed formulae for heavy quark decays. Effects of the final quark mass are taken into 
account. Explicit expressions are obtained for the t — > b and b — > c transitions for both the 
frozen coupling approximation and in the QCD running coupling case. 



1 Introduction 



A common feature to many processes in QCD is the presence, in the perturbative expansion, 
of large double (Sudakov-like) logarithms at the threshold. Resummation of large infrared loga- 
rithms in form factors and shape variables is essential in order to predict accurate cross sections 
in many phenomenologically relevant processes (see, f.i. [1-5]). In semi-leptonic heavy quark 
decays — > qf I u, threshold regions are characterized by the presence of two different scales 
mx <C Ex, where mx and Ex are the final hadron invariant mass and energy, respectively, 
originated by the final quark qf. The perturbative expansion is spoiled by logarithms of the 
ratio of the two scales. Those need, therefore, to be resummed. 

Such logarithms are organized as a series of the form [6, 7]: 

°° 2n O 2 O 2 O 2 

VV Crifc <(Q)log fe ^V= c 12 a s (Q) log 2 \- + Cll a s (Q) log \- + (1) 

n=l k=l A A A 

+ c 24 a 2 s (Q) log 4 ^ + c 23 a 2 s (Q) log 3 ^ + ■ ■ ■ , 
m x m x 

where a s is the QCD coupling constant and Q is the hard scale Q = 2E X - The leading term 
is the double logarithm a s (Q) log 2 Or. A similar double logarithmic structure is present in 
many other processes like deep inelastic scattering (DIS), heavy quark fragmentation, Drell- 
Yan annihilation, Higgs production, and so on, the argument of the logarithms differing per 
observable and per process. 

A universal resummation formula valid at all perturbative orders for a decay of a heavy quark 
into a massive quark qf plus a non-hadronic state, with a final state jet-like structure, has been 
recently obtained [8]. 

In QCD resummed formulas, the running coupling is integrated over all gluon radiative 
momenta from the hard scale down to zero, hitting the Landau pole. A prescription has to be 
assigned to give a meaning to the formal resummed expressions. 

One possible solution is the use of an additional prescription for the contour integration 
in A^-space, in the inverse Mellin transform from A-space to x-space, the so-called minimal- 
prescription (MP) [9]. This prescription provides a formula which is the asymptotic limit of the 
expansion, furthermore it is renormalon free and the truncation of the series at the minimum 
term originates an exponentially suppressed difference between the truncated expansion and the 
full MP formula. 

The aim of this work is to analyze the perturbative resummed distributions in the parton 
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subprocess for both massless and massive final heavy quarks and to explore the feasibility of the 
MP regularization scheme. 

We apply the resummation formulas to the case of t — > b and b — > c, as a working example 
to implement this regulation method. In literature, the b — > c inclusive semileptonic decays are 
widely discussed, also in the contest of effective theories, in order to improve comparison with 
the newest data (for a review see f.i [10] and refs within). The t — > b case has been discussed in 
QCD resummed formulas [11], with different dynamical variables. 

The assessment of a perturbative reliable and singularity safe form factor is the first step 
towards a sound phenomenological approach and it is also needed for comparison with QCD 
based effective theories like, for example, SCET. 

The papers is organized as follows: in section 2 we recall the resumming formulas for the 
final massless and massive quark in the Mellin space, and recast them in a more transparent 
notation; in section 3 we move to the physical space and analyze the feasibility of the MP for 
the massive case. In section 4 we study the frozen coupling approximation, while final plots with 
the QCD running coupling and conclusions are presented in sections 5 and 6. 

2 Threshold resummed Jet Distribution in Mellin space 
2.1 Massless final quark 

Before considering the case of a massive final case, let us recall the expressions and the variables 
for the resummed jet distributions in the massless final state [5, 12-15]. 

Let us consider the decay driven at a partonic level by an heavy quark decaying into an 
approximately massless final quark, plus non-hadronic states, as, for instance, the decay B — > 
X u lu or the radiative decay B — > X s ^. Threshold resummation is typically performed in Mellin 
space; the threshold limit corresponds to iV — > oo and threshold logarithms a™ log™ N can be 
factorized into a form factor Jn, which has the exponential form: 



Jn(Q 2 ) is the massless jet distribution, that gives the probability that a massless parton produced 
in a hard process with a hard scale Q fragments into a hadronic jet of mass mx 



Jn(Q 2 ) 



e /^(Q 2 ) ) 



(2) 



m\ = (l-x)Q 2 



(3) 
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The Mellin or iV— transform is defined as 

J N {Q 2 ) = I' dxx N - 1 J(x;Q 2 ) (4) 
Jo 

In the limit a s — > 0, the mass distribution reduces to a spike corresponding to the (zero) 
parton mass. In the limit x — > 1 we drift away from the perturbative regime. If x — 1, the 
truncated perturbative expansion becomes unreliable. It is possible, however, to be able to use 
a perturbative resummed expression at all orders in the Mellin space which reads [5, 12-15]: 

Jn(Q 2 ) = exp f [x N -> - 1] ( f Q2{l ~ X) d ^A [ a , (k 2 ± )} + B [a s (Q 2 (l - x))] + 

Jo 1 _ x y Jq 2 {i-x) 2 fc ± 

+D[a s (Q 2 (l-x) 2 )} | (5) 

The functions A (a s ), B (a s ) and D (a s ) have a perturbative expansion: 

A(a s ) = A 1 a s + A 2 a 2 + --- , B (a s ) = B 1 a s + B 2 a 2 s + ■ ■ ■ , D {a s ) = D x a s + D 2 a 2 s + ■ ■ ■ . (6) 
The known values of the coefficients A i: B; L and Di are given in [14, 16, 17]. 

A (a s ) describes the emission of partons which are both soft and collinear, B (a s ) describes 
hard and collinear partons while D (a s ) partons which are emitted soft at large angles. A(a s ) and 
B(a s ) are related to small-angle emission only. They, therefore, represent intra-jet properties[16, 
17], while the function D(a s ), being related to soft emissions at large angles, is a process- 
dependent inter-jet quantity. 

While the validity of the resummed formula goes beyond our case of semi-leptonic heavy quark 
decays, holding for accounting of threshold logarithms in several other processes, the specific 
structure of (5) can vary, depending on the specific process and on the particular observable under 
exam. For instance, there are corresponding results for the DIS structure functions -Fi j2 ,3(a;, Q 2 ), 
where Q 2 represents the resolution scale, or for the Drell-Yan cross section da/dQ 2 , where Q 2 
stands for the invariant mass squared of the lepton pair [13, 18]. 

In order to illustrate how to interpret the universality of formula (5), let us consider the order 
a s decay t — >■ bW g (where W and g are a real W boson and a gluon, respectively) and examine 
the distribution in the energy of the final 6-quark, that is in the variable Xb = 2Eb/mt. Once 
considering the distribution in Xb, rather than in x as in the present paper (see definition (3)), 
we are dealing with a different observable-and therefore a different kinematical parametrization 
of the threshold region. Now the threshold region is reached when Xb — > 1, a limiting point 
where there is no gluon emission to change the light quark energy. We expect only the emission 
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of soft gluons and we do not need any B(a s ) contribution in (5), since this function contains 
collinear radiation associated with the light quark. Formula (5) still holds, but without the 
B(a s ) termfll]. In the present work, collinear gluons described by B(a s ) are allowed, since the 
energy change of the light quark, still in the jet after the gluon emission, does not affect the 
distributions in x, related to the invariant mass of the jet. 



The exponent of Eq. (2) can be expanded in a function series of the form [13]: 

oo 

f N (Q 2 ) = log J N (Q 2 ) — Lgi (A) + $>?<7 n+2 (A) = L 9l (A) + g 2 (A) + a s g 3 (A) + 
where 



n=0 



and (3 



A = A, «s(Q 2 ) L, 
(11/3 N c - 2/3 n F )/{Au). 



L = logiV 



(7) 



(8) 



The first exponential term Lg\ (A) = L Xl^Li 9i,n^ n resums the leading logarithms (LL); by 
adding the term g 2 (A) = Yl^=i fi^nA™, also next-to-leading order terms (NLL) are taken into 
account and resummed, and so on. 



The functions gi (A) have a power-series expansion: 

oo 



(9) 



n=l 



They are all homogeneous functions: g«(0) = 0. This property insures the normalization of the 
form factor Jn=i = 1- The functions g± and g 2 become singular, signaling non-perturbative 



effects, at A = 1/2, that is at iV = exp[l/2/3 a: s (/! 2 ) 
[19,20]. 



/x 2 /A 2 . Explicit expressions are given in 



The leading logarithmic term in Eq. (5) is 
J N ~ exp [L g x (A)] ~ exp 



a s L 

2 



exp 



(10) 



Function series like Eq. (7) appear in other processes as well as, for instance, DIS and Drell-Yan 
with the same towers of threshold logarithms. Additional terms, however, due to soft-gluon 
radiation collinear to the light initial-state parton, in the DIS case, and to two light initial-state 
partons, in the Drell-Yan case, do slightly modify the form of g±, giving, for the leading term of 
the resummed quark coefficient functions C^ IS and C^y [18] respectively : 



C£ /s ocexp[Ls? /s (A)] ^exp 



C N DY ocexp [Lg? Y (\)} ~ exp 



2^1og 2 iV 

7T 
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We plot in fig. 1 the form factor in the massless case, for the £—>■£> and the b — > decays, 
in order to show the different slopes in the two cases, that will be maintained in the massive 
case, affecting the regularization procedure. We have set the scale /i in equal to Q 2 . The 
continuous and dashed lines represent the NLL and NNLL contributions respectively. The 
strong dependence on a s values is shown: light and thick lines are referred to different values of 
a s . The NNLL curve stands below the NLL one, due to the inclusion of g$, that is a negative 
decreasing function within the considered range. We also observe a relatively strong dependence 
on a s . As said before, the effects of the Landau pole start appearing at A = 1/2, that is towards 
N ~ 10 6 for the case of the top decay, and for N ~ 10 2 for the case of bottom decay; there, the 
Mellin form factor starts to oscillate and the expressions are no longer predictive. 




J N 



500 1000 1500 10 20 30 40 



Figure 1: Form factor in N space for massless final quark. Left figure: the t — >■ b case: a s = 0.11 
(light lines), a s = 0.12 (thick lines). Right figure: the b — > c case, a s = 0.20 (light lines), 
a s = 0.22 (thick lines). In both figures the continuous lines represent NLL contributions, the 
dashed lines NNLL contributions. 

2.2 Mass-corrected Jet Distribution 

Let us briefly summarize the results obtained for the case of a massive final quark qf [8] . The more 
massive is the radiating qf, the less radiation has to be emitted in the decay; as a consequence, 
the typical Sudakov effect, namely the suppression of non-radiative channels and the broadening 
of sharp structures, are expected to be less pronounced for the massive channels. In principle 
one has single-logarithmic corrections, which are not strong enough to shift the peak of tree-level 
distributions. 

In Ref.[8] it was demonstrated at NNL order (and conjectured to be also valid at accuracy 
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beyond NNL) that the massive jet function can be factorized in momentum space, as: 

Jn(Q 2 ; r) = J N (Q 2 ) 5 N (Q 2 ; r) (12) 
Jn(Q 2 ) is the massless jet distribution and Sn(Q 2 ; r) is the mass-correction factor which reads: 

rl T r(7V-l)_i f r m 2 (l-x) jj.2 

5 N (Q 2 ; r) = exp / dx—^ -\ - / [a s (k 2 ± )} - B [a s (m 2 (l - x))} + 

J0 l - X y Jm 2 (l-x) 2 fc _L 

+ D[a s {m 2 (l-x) 2 )} I, (13) 



x and r are defined as 



2 2 2 



where Q is the hard scale of the process and m is the mass of the emitting quark. We assume 
the quark mass to be much smaller than the hard scale, in order to have fast-moving charges and 
to preserve a jet structure. We indicate both the mass corrected and the massless jet functions 
with Jjv; they are distinguishable since the the massive one bears a dependence on r. 

Eq. (13) has a simple physical interpretation. The parameter ^ — 1 is multiplied by r on the 
r.h.s. of Eq. (13), implying that mass effects become "visible" only for large 

N > - > 1. (15) 
r 

In this case, there is enough resolution to "see" the quark mass, which tends to suppress the 
collinear effects, related to the A and B terms. At the same time, soft radiation not collinearly 
enhanced, described by the function D and characteristic of massive partons, does appear. Let 
us also note that, since the jet mass is an infrared (i.e. soft and collinear) safe quantity, 5n = 1 
for r = 0. In the limit r — > 0, the well-known massless result is recovered. 

The mass-correction factor has the same structure than the massless case [13] 

5 N (Q 2 ;r) = e M<?V) ; (16) 

where the exponent has a double expansion of the form: 

oo n+1 

F N {Q 2 ;r) = 9 (N — 1/r) £ £ F nk a n s log fc (iVr), (17) 

n=l k=l 

with F n k numerical coefficients. The exponent can be expanded in towers of logarithms as: 



F N (Q 2 ;r) = L d, (p) + ]T a?d n+2 (p) 



n=0 
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= Ldi (p) + d 2 (p) + a s d 3 (p) + a]d 4 (p) + • • • , 



(18) 



where 

p = f3 a s (p 2 ) L, and L = 9 (N - 1/r) log (N r) . (19) 

The scale p = 0(m) is a renormalization scale of the order of the quark mass m. The over- 
all factor 9 (N — 1/r) comes from the step approximation of the moment kernel and avoids 
modifications for small iV of the massless behavior, in agreement with the physical intuition. 
Furthermore it ensures the correct massless behavior in the r — > limit. Analytic continuation 
to the complex iV-plane can be made by omitting such factor and fixing the correct interval in 
physical space. 

By truncating the above series expansion, one obtains a fixed-logarithmic approximation to the 
form factor 5n- Functions di(p), which represent the mass effects, can be obtained from the 
standard ones gi(\) of the massless case [7] by means of the replacements: 

2 2 

A(a s ) ->■ -A{a s ); B{a s ) -> - B(a s ); D{a s ) -> D(a s ); log ^ -> log^; A ->■ p. (20) 

It is worth observing that mass effects induce a similar structure to the massless one, involving 
changes of sign of the collinear functions A and B, with the rescaling Q — >■ m. The explicit 
expressions for the functions di are listed in [8]. 

Let us now examine the behaviour of the jet function as given by the equation (12). Through- 
out the paper we fix the hard scale of the process Q to the mass of the decaying quark, that 
is to m t in the case of top decays, and to rrib for b decays. The correction factor 5n{Q 2 ', r) 
is a function increasing with iV; in order to produce Jn{Q 2 ', r ) h has to be multiplied by the 
massless form factor Jn(Q 2 ), at values of iV > 1/r. In the case of the top quark decay, the 
increase is very slow; it starts at 1/r ~ 2xl0 3 and the distribution only doubles when iV reaches 
~ 6x 10 4 , continuing slowly, until a fast increase before values of ~ 2x 10 m 2 /A 2 ~ 10 6 , where it 
reaches the peak 1 . That implies, as expected, that mass addition does not modify substantially 
the massless distribution until very large values of N. The left hand plot in fig. 1 practically 
coincides with JN(Q 2 , r ) defined in (12), at the same value of a s . 

On the contrary, in the b — > c case, the increase starts earlier (1/r ~ 2x 10) and it is much faster 
(it doubles one at iV ~ 10 2 ), reaching the fast increase and the peak around iV ~ m 2 / 'A 2 ; effects 
are much more sizable. NNLL corrected curves grow faster than the NLL curves, as shown in 
fig. 2. 

1 Differences between the mass correction factor at NNLL and NLL order are that the NNLL corrected one 
peaks much faster and at an earlier point in TV; the increase, however, still occurs at the same order of magnitude 
and it does not introduce substantial changes on the distribution. 
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Figure 2: Form factor in N space for b — > c at a s = 0.219: the continuous and dashed lines 
represent NLL and NNLL contributions, respectively. Light and thick lines refer here to massive 
and massless final quark, respectively. 



3 Threshold resummed Jet Distribution in physical space 

Even if the iV-moment expressions of the jet function are physical quantities, their measurement, 
especially for large N, is difficult. It is therefore convenient to perform the inverse Mellin 
transform back to momentum space. Given the Mellin transform Jn defined as in (4), its inverse 
transform is 

i rC+ioo 

J(x;Q 2 ) = M'^Jn-x] = — / dNx~ N J N (Q 2 ) (21) 

2m Jc-ioc 

The inverse transform of the product of two generic /jv and g n is the convolution of the two 
inverse functions f(x) and g(x): 

M- 1 \f„g N ;x]= f f(-)g(u)- (22) 
J x \uJ u 

3.1 Massless final quark 

There are two possible ways of obtaining the J(x; Q 2 ) distribution from the inverse Mellin 
trasform. Each of them has its own peculiarities, since we are dealing with truncated expressions. 
We have used both in order to compare the results and increase their reliability. 

One way is to use an analytical expression for the inverse Mellin transform (21). Indeed, the 
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massless form factor Jn(Q 2 ) is 



J(x-Q 2 ) = - x -^-|0(l - x -e) E (x; Q 2 ) } at e ->■ 0, 



(23) 



where £(x; Q 2 ) is the inverse Mellin transform of Jn/N. 



The 0(1 — x) function ensures the unitary normalization of the distribution in the interval 
(0, 1) and it can be omitted in the massless case since the function is regular at the origin. 

We have, at NNLL [19] that : 



where 



1 + F»' L ^(1 - F? L ) + - F 2 (l) (^(1 - F? L ) - ^(1 - F» L )) 



r(i - L ) 

e l g 1 (l3 a s l)+g2(fSoa s l)+a s g 3 (j3 a s l) 

r(l - gi((3Qa s l) - p a s lg[(Poa s l)) 



[1 + /3 o a s fi'2(/ 3 oa s 0V ; (l -0i(A)a s O - A)"* ^(A)"^)) 



+-F 2 (Z) (^ 2 (1 -<?i(/W) -/W<7i(/W)) - - <?i(/W) -/W<7i(/W))) 



(24) 



*i>(0 
Ff 2L (/) 

^2(0 



lgi(f3 a s l) + giifootsl) + a s g 3 ((3 a s l), 
gi(Poa s l) + f3 a s l g[(f3 a s l) 
(3 a s g' 2 (P a s l). 

2/?o«s ffi (/W) + /3 2 a 2 / s"(/W) . 



Here T is the Euler Gamma function, ip(x) = d log T(x)/dx, the digamma function, and 
I = — ln(— lnx). Note that / — >■ — ln(l — x) when x — > 1. 

Expression (24) can be rewritten in a synthetic way by evidencing the NLL part, that is as 

p l 9i(t) + 92(r) 

^• < ^ = W^m] SE (25) 

with 

r = /W, h 1 (r) = ^(Tg 1 (r)) (26) 

and 

= K x e °-*>M jl + o a 9 g' 2 (T)1> [1 - /n(r)] + ^A> « s /ii(r){^ 2 [1 - /ii(r)]-^' [1 - ^i(r)] } 

(27) 



9 



Here, K 1 is a normalization factor such that <5£ — > 1 when / — > (or rr — > 0). 



Another possibility is to obtain the inverse Mellin transform numerically, by integrating 
Eq. 5) at next-to-leading order. 

This numerical integral is not straightforward, since, as we have seen, the gi are singular in 
A and their singularity reflects into N. In other terms, the numerical distribution is not real for 
any value of N because of the integration over the Landau pole. An exact numerical evaluation of 
the inverse transform then requires a prescription for the pole. We use the Minimal Prescription 
(MP), on a suitable path to the left of all the singularities [9]. We have compared the analytical 
distribution (23) with the distribution obtained numerically. The two curves show a very good 
agreement, although they differ slightly around x ~ 1, since the analytical ones reach the peak 
and start oscillating earlier- 



3.2 Mass correction factor 

In analogy to the massless case, the mass-correction factor in physical space is obtained by means 
of the derivative of the inverse Mellin transform of Sn/N: 

^ ;0 W) = -4{£J^-MQW)}, (28) 

where c is a (real) constant chosen in such a way that the integration contour lies to the right 
of all the singularities of 5n- By defining 

5 Nr = 5 N (29) 
and by changing variable from N to v = Nr, we obtain: 

% ;QW) - -*± {H^^W™ 1 )}. (30) 

After this change of variable we can neglect the 9[u — 1] and make the analytic continuation 
in the complex N-space. 

We can therefore use the results in [19] to obtain the correction factor in physical space in 
NNLL approximation: 

5 (x; Q\ m 2 ) = - x -jL{e(l- x - e )A (x; Q 2 , m 2 ) } at e ->■ 0, (31) 

where: 

J'd 1 (T')+d2(T') 

A ^^')' r[i- M i a (32 > 
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is the resummed partially integrated (or cumulative ) form factor and the 8(1 — x — e) ensures the 
unitary normalization of the distribution in the interval (0, 1). As already observed, this term 
can be omitted in the massless case since the function is regular at the boundary, but protects 
the mass correction factor which is not a regular physical distribution. 

In Eq. (32) we have defined 

l' = - log(- logx 1 ^) t' ee /W- (33) 

and 

h^r') ee A [Xd^r')] = rfx(r') + Xd' 1 (r'). (34) 
SA is a NNLL correction factor which can be set equal to one in NLL: 

SA NLL = 1. (35) 

Its NNLL expression reads: 

SA = (36) 

with 

S = e a ^') jx + p a s d' 2 (T')i>[l - hir')} + ±/? a s h[(r'){^ [1 - MOW [1 - fciM] } J- 

(37) 

IYx) is the Euler Gamma function and 



^(x) ee Alo g r(a;) (38) 
ax 



is the digamma function. 



It is convenient to approximate the argument of the inverse Mellin transform for y = 1—x <C r 
by the expansion: 



Hi=i-f+o(g, 



(39) 



so that 



d f f cr+lo ° dv 



i (y, Q 2 ,m?) = (1 - „) £ | jf ^ [ 1 - l\ K (Q\ »') j . (40) 

Note that the r.h.s. is positive only for y < r, implying that the linearization above shrinks 
the domain of y from the unitary interval (0, 1) to the much smaller interval (0, r). The correction 
factor in physical space is therefore the inverse Mellin transform of S u with respect to (1 — -) 



— V 
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In this case it is useful to employ the limit definition as for the plus-distribution defined in 
Ref. [7], such as: 

d 

dy 



5(y;Q 2 ,m 2 ) = lim (1 - y) A [o (y - e) A (y; Q 2 , m 2 ) } 



Finally: 



t' = f3 a s L 



and 



1- V - 



L — — log - log 

A further approximation step can be made in order to obtain the final result: 

,2~ 



log 



l- V - 



y - + o( y -). 

r \r z 



Finally the resummed expression in physical space reads: 

5 (y; Q 2 , m 2 ) = lim + (1 - y) |- [o (y - e) A (y; Q 2 , m 2 ) } 
where A (y; Q 2 ,m 2 ) is given by Eq. (32) and 



L — 9(r — y) log 



(41) 



(42) 
(43) 

(44) 



(45) 



(46) 



We have limited the domain to y < r with a ^-function 2 . 

We are now ready to perform the convolution in order to obtain the physical distribution. 
The physical form distribution is obtained by the Mellin transform of Eq. (12), that is by 

rc+ioo dN 



J 



/■C+lOO TAT /-C+iOO TAT 

(x; Q\ r)= ^ x~ N J N (Q 2 ; r) = %L X ~ N J N {Q 2 ) 5 N (Q 2 ; r) (47) 



This integral is not straightforward, since, as we have seen, the Qi are singular in A and their 
singularity reflects in N. 

J (y; Q 2 , r) can also be computed analytically, by the convolution of the inverse Mellin trans- 
forms of Jn(Q 2 ) and Sn(Q 2 ; r). 



J(x;Q 2 ,r) = l 1 ^J(z;Q 2 )5{^ z ,Q 2 ,r) 



(48) 



2 As suggested in Ref. [8] a smooth approximation to the Theta <8> Log function form can be given by: 
L ~ — log [l — (1 — y) 1/V ] . In fact these functions agree at the first order approximation. 
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where J(y; Q 2 ) and S (y; Q 2 ,r) are given by Eq. (23) and Eq. (31), respectively. Therefore, we 
obtain the following analytical expression: 

J(x;Q 2 ,r) = f d ^J{z-Q 2 )5(*-Q\r) (49) 

J X Z ^ Z ' 

= hm /^z- {-5(l-^-e)S'(- ; g 2 ) A(z;Q 2 ,r) + 6(1 - z - e)^ (-; Q 2 ) A' (z; Q 2 ,r)\ 

e^0+ J x Z < Z Z J 

where we have the Dirac delta function 5(1 — z — e) = — d6(l — z — e)/dz. Let us observe that 
Z(y,Q 2 ) ->■ when y ->■ 0. 



4 Frozen Coupling approximation 

The frozen coupling approximation means neglecting the variation of a s with the scale. We 
first look for solution in the frozen coupling approximation; in the massless case the resumming 
formula at NNLL gives: 

f 1 x 1 ^' 1 - l r r Q2 ^ x) dk 2 

log J N (Q 2 ) = / dx — \ / — [Akx, + A 2 a 2 + A 3 a 3 s + ...] + 

JO l — X l jQ2(l- x )2 K 

+ B ± a s + B 2 a 2 + ... + D iag + D 2 a 2 s + ... + } 

f 1 x N ~ l - l r l 

~ / dx— <^ (A^s + A 2 a 2 s + A 3 a s s ) In + 

Jq 1 — x < 1 — x 

+ (B 1 + D 1 )a s + (B 2 + D 2 )a 2 } 

(50) 

In the frozen coupling approximation, fl — >■ 0. After integration in we expand in A = 
/3o« s £ where L = log TV. At the lowest order, in the massless case, we have: 

91 = "^-A (51) 

/ B x D x A llE \ 

92 = h^-^-^o-r (52) 

93 = h^7^— r (53) 

We can easily find the corresponding di for the mass correction term by the substitution 
Ai -» -A i: Bi ->■ -Bi and A ->■ p. 
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Disregarding the NLL terms we have in the frozen coupling limit: 

log J N (Q 2 ,m 2 ) = f N (Q 2 ) + F N (Q 2 ,m 2 ) 
~ Lg 1 + L r d\ 

= -^(XL-pL r ) (54) 

where p = a s f3oL r with L r = log Nr. Then the leading behavior for the mass corrected formula 
is determined by: 

log J N (Q 2 ,m 2 ) ~ A 1 a s logiV logr (55) 

The divergent double logarithmic behavior for iV — > oo of the massless case is replaced by a 
single logarithm times a new regularizing term logr. This term is negative and restores a finite 
limit when x — > 1 in the physical space. The peculiarity of applying the MP is that due to the 
milder singularity, we no longer have the factorially growing spurious contributions described in 
[9] generated by neglecting certain subleading terms when the moment space formula is turned 
to an rr-space formula. The resummed massive case distribution is then a regular function in 
the limit x — > 1. We have therefore found that, also in the massive case, the resummed formula 
is void of unwanted spurious ambiguities. 

i dr 

r i i dr 

r dx R(x)= | dz 

5000 r 

4000 : 
3000 - 
2000 : 
1000 - 



0.990 0.992 0.994 0.996 0.998 1.000 * "o.80 0.85 0.90 0.95 1.00 

Figure 3: Frozen coupling approx.: top decay jet rates (on the left) and partially integrated jet 
rates (on the right). Comparison between massless (continuous line) and massive (dashed line) 
distributions. 

It has been already mentioned that there are two ways, analytical and numerical, to compute 
the inverse Mellin transforms of Jn{Q 2 ', r). By numerical method we mean the direct numerical 
integration in Eq. (47); by analytical method we mean to use the approximated analytical 
expression for the convolution in Eq. (50). In the frozen coupling case, the and di are linear in 
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Figure 4: Frozen coupling approx.: bottom decay jet rates (on the left) and partially integrated 
jet rates (on the right). Comparison between massless (continuous line) and massive (dashed 
line) distributions. 

A and therefore the numerical path does not include the Landau pole; the numerical integration 
becomes therefore exact. 

In figs. 3-4 we compare the (normalized) resummed massless and massive jet rates, in the 
frozen coupling approximation, for top to b and for b to c decays, respectively. We obtain the 
same results by calculating the jet factor with both numerical and analytical methods. 



5 Resumming with a running coupling constant 

By releasing the frozen coupling approximation, we can calculate the distribution with a running 
coupling. We need in this case a regularization procedure. Indeed, in the massive case, we have 
four poles on the real axis of the Mellin complex plane respectively for A = |, A = 1 and for 
p — |, p — 1. These poles arise from the two logarithmic structures of the massless and massive 
correction formulas. The important point is that the massive poles stand to the right of the 
massless ones, therefore allowing to use the MP procedure of Ref. [9]. We have integrated 
numerically Eq. (47), over the path made by two straight lines parallel to the negative real 
axis, closed by a half-circle centered around the origin and crossing the positive axis between 
the origin and the first Landau pole. The integrated expression is computed at NNLL order, 
releasing the frozen coupling approximation; the coupling runs over the whole integration range. 

In fig. 5 we compare the resummed massless and massive jet rates for the case of t to b decays. 
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i dr 

r dx 

500 r 




— i I i 1 1 1 1 1 1 1 1 1 1 1 — \ 1 1 1 1 1 1 — « — i 

0.9995 0.9996 0.9997 0.9998 0.9999 



Figure 5: Top decay jet rates: the massive case (dashed line) compared to the massless case 
(continuous line) at NNLL; a s is set to the top mass scale, a s = 0.11. 

All the plots in this section are not normalized. The massless plots coincide, in the considered 
range, with the ones obtained by the analytical distribution of Eq. (23). Differences between 
massless and massive case start, as expected, approximately for values of x ~ 1 — r ~ 0.999 3 . We 
have only listed the NNLL plots, since they do not differ, substantially, with respect to NLL ones. 
The MP does not avoid approaching to the essentially non-perturbative regime at x ~ 1, where 
the plots start to oscillate and, therefore, a physically motivated treatment of non-perturbative 
effects has to be introduced. We find that this physical non perturbative cut-off can be put, in 
practice, equal to 1 in top decay. 

In the b decay case (fig. 6), the two curves start differentiating at x ~ 1 — r ~ 0.9, and 
the difference is visible, with the massive case being less divergent, as expected. The rates in 
the curves are not normalized, and the NNLL plot presents relevant difference from the NLL 
one; the addition of NNLL terms renders the plot closer in the rising behavior to the massless 
plot. At NNLL order we deal not only with double logarithms, but also with single logarithms 
that shift the position of the minimum 4 . The effect of the running coupling, affecting in a very 
distinct manner the massless and the massive case, can be observed by comparing with fig. 4. 
In order to facilitate the comparison with the frozen coupling case, we report both the massless 
and massive distributions, not normalized, with frozen and running couplings, in fig. 7. 

3 We have, approximately, z N ~ x — 1 ~ —8(1 — z — 1/N) [13], and we can set x < 1 — 1/N ~ 1 — r. 
4 scc f.i. formula (62) in [7]. 
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Figure 6: B decay jet rates: the massive case (dashed line) compared to the massless case 
(continuous line) at NNLL; a s is set to the beauty mass scale, a s = 0.219. 



6 Conclusions 



In this paper we have considered the application of the MP to the resummed jet function dis- 
tributions for heavy quark decays in massless and massive final quarks. We have used a QCD 
resummation formula in a s , that takes into account contribution from large collinear and soft 
logarithms near the threshold region [8]. Benefits of large threshold logarithms resummation 
are restoring predictive power to the perturbation series and increase theoretical accuracy, f.i. 
by a reduction of scale uncertainty. Accuracy is obviously increased by including higher order 
terms in the exponent. We have considered NNLL order corrections. In the case of massive final 
quarks, we have analyzed the applicability of the Minimal Prescription scheme to deal with the 
unavoidable problem, connected with any physical application of QCD resumming formulas, of 
the integration over the Landau pole. 

We calculated the form factors in Mellin and physical space, in order to extract the main and 
universal features of the perturbative distributions, common to all processes at the threshold. 
We consider this a necessary first step to approach the phenomenological study of specific decay 
rates in the QCD resummed framework. In order to compare with data, the further step is 
to match the distribution with a fixed order coefficient function and include non perturbative 
effects [21] 

As a final remark, let us notice that a very interesting application is to the inclusive semi-leptonic 
b decay into c. In such partonic three body decay, one can combine the resummed with the full 
triple differential distribution. 
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0.998 



1.000 



Figure 7: Comparison between the frozen coupling (lighter/red lines) and the running coupling 
case (thicker/black lines) in B decay jet rates. The dashed (continuous) lines refers the massive 
(massless) case, at NNLL and with a s = 0.23. 

We have considered top quark decays, for the recently increased interest due to the large amount 
of data available at the hadronic machines. One possible application of the present results is to 
explore the effects of the resummation on the calculation of helicity fractions of the W boson 
from top quark decays. These are presently measured with increasing precision by both the CDF 
and the DO collaboration at the Tevatron. 
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8 Appendix 

In this appendix we enlist notations and actual values used in the paper. 

The functions gx and #2 introduced in section (2.1) have the following expressions [20]: 



91 ( A; W 



) 



2/3 A 



A x 1 



[(1 - 2A) log (1 - 2A) - 2 (1 - A) log (1 - A)] ; 
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= [log(l - 2A) - 2 log(l - A)] + ^ [log(l - 2A) - log(l - A)] + 

Po 

- [log 2 (l-2A)-21og 2 (l-A) + 21og(l-2A)-41og(l-A)] + 



2ft 



*ft 



A 1 



log(l - 2A) + ^ log(l - A) + [log (1 - 2A) - 2 log (1 - A)] log ^ 



2/?c 



The NNLO function g 3 [19] reads: 



93 ( A; b 



A, 



A 



A 



2ft 

^2 

"4/3 3 



1 - 2A 1 - A 
2A 2A 



A^2 



2/3, 



4A 



<2 2 



A 



+ 



2ft 
Ax/3 2 



1-2A 1-A 

2A 2A 
I _^--_ + 2tog(l-2A)-41og(l-A) 

log(l-2A) 2 log (1-A) 3A 3A 
— — .-• l - A + 1 - 2A ~~ 1 - A_ 



2ft 



1 - 2A 

2 



l log 2 (l-2A) _ log 2 (1 - A) + log(l-2A) 



+ 
+ 



1 - 2A 



DiPi pog(l-2A) 2A 1 f log (1-A) J_ 



1-A 
A 



1 -2A 



2ft 

D 2 A 



1 -2A 



1 -2A_ 

£ 2 A Ag% 
?o 1 - 2A ftl-A 2 



ft 
4A 



1-A 
A 



" ft 

A 2 1e 



,, . .. _ L l -2A 1-A 
log(l-2A) log (1-A) | 1 



+ 
1 



A) 
A 1 



1 -2A 
1 1 



2/?o 
A) 



1 -2A ~ 1 - A 
2A 2 A 2 
1-2A 1-A 
2A A 



1-A ' 1-2A 1-A 
Dg E 2X B llE X 



1-A 



+ 



+ 



+ 



1 -2A 

.2 



log 2 

g Q 2 ft 



'0 



1-A 

A 



+ 



A 



1-2A 1-A 



1 - 2A 1-A 



^~ -Pi A /z~ -pi a , fJr 

, .. .. , o- ;„ i - 2A og g 2 ; (l 1 a og g 2 + 

Alog(l-2A) A log (1-A) | A 



Ex A 



1 -2A 



1-A 



1-2A 



i V 



ft 

- r ^ + ^log(l-2A)-log(l-A) 

Arbitrary constants have been added to the function g 3 in order to make it homogenous. The 
quantity 7^ = 0.577216 ... is the Euler constant and ( (n) is the Riemann zeta function, 
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C (2) = 7r 2 /6 = 1.64493. The functions g 2 and g 3 depend on the renormalization scale /i, while 
gi does not. 

The known values for the resummation constants defined in section (2.1) read: 



A 1 
A 2 



7T 
TT 2 

Cf 

7T 3 



245 



£(2) 
2 

11 



5 

67 



/209 



/55 z(3) 



+ 



108 



S 2 



3CV 



732 



4" 


7T 


7T 2 


C.4 


- c F 








V32 
'247 




^432 


_C F 


1 

7T 


c F 


C A 


TT 2 





3155 



11^(2) + ^(3)1 f 



864 12 



+ -z(3)--z(2)) - 



z(2) 



55 
108 



9 

T 



z(2) 



54 



where = N c = 3 is the Casimir of the adjoint representation. 

The knowledge of the quantities A±, A 2 , B\ and D\ is needed for resummation at next-to- 
leading order. 

The coefficients A\, B\ and D\ are renormalization-scheme independent, as they can be 
obtained from tree-level amplitudes with one-gluon emission. The higher-order coefficients are 
instead renormalization-scheme dependent and are given in the MS scheme for the coupling 
constant 5 . 



5 A discussion about the scheme dependence of the higher order coefficients A2, £>2, etc. on the coupling 
constant can be found in [14]. 
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The coefficients dj, denned in section (2.2), are : 



di(p) 
d 2 (p) 



, (l-2p)log(l-2p)-2(l-p)log(l-p) 

2 Pop L 

^-log(l-2p)-|ilog(l-p)--^ 

z Po Po ^ Po 



+ 



4 A) 3 



log(l-2p)-2 1og(l-p) 

+ 



+ 



2 log(l - 2p) + log 2 (l - 2p) - 4 log(l - p) - 21og 2 (l - p) 



log(l-2p)-log(l-p) 



2 A) 



p/ 



log(l-2p)-2 1og(l-p) log-^. 



For the NNLO function <i 3 we obtain: 
D 2 P 



+ 



+ 



+ 
+ 



+ 
+ 



+ 



A) l-2p 
5 2 p 



2p £>i A 



2p + log(l-2p) 



A 1-P 

5! A 



+ 5x7^ 



P 



+ 



P 



2 /V 
2 

Ai A 

4pV 
Ai gi 7ff 

A 2 
A 1 /3 1 2 

2 A) 4 

,2/ 



+ 



1-P 

log(l - p) 



1-p 1-p 

2p p 

l-2p 1-p 
3p 



+ 



As 



2/V 



l-2p l-2p 



+ 



+ 



l-2p 1-p 



+ 



l-2p 1-p 
4p p 



3p + log(l-2p) 2 1og(l-p) 



l-2p 1-p 
2p 2p 



+ 



1 -2p 

A 1 7T 2 



1-p 



12 



4p 



l-2p 1-p 



+ 
+ 



1 -2p 1-p 



+ 2 1og(l-2p)-4 1og(l-p) 



2p 



P , log(l - 2p) log(l - p) 



1 -2p 1-p 1 -2p 1-p 

^-log(l-2p) + 10g(1 - 2p) 

l-2p 1-p &K l! l-2p 



+ 
+ 



log 2 (l - 2p) + 2 _ 2 log(l - p) log 2 (l - p) 



2(l-2p) ' -" - - 1-p 
D\ p p 2 B x p p? A 2 



1 -h — 

A, l-2p ~~°m 2 ' /3 1-p ° g ^2 ^2 



AiTg 

A 
Aip\ 

A 3 



2p i> 



1-p 

P 



+ 



P 



l-2p 1-p 



, P 



l-2p 1-p 
P P 



i P 



1 - 2p 1-p 



| log(l-2p) | plog(l-2p) | 



l-2p 



log(l - p) 



P log(l - P) 

1-p 



, p 

log^ + 
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+ 



A 

2 A) 



2p 2 



l-2p 1-p 



, 2 A* 
log — ». 

rrr 



The coefficients of the QCD /3-function in our normalization have been given in [7] . 
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